y=dsolve('x*D2y+(1-n)*Dy+y==0','y(0)==0','Dy(0)==0','n')